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A Mathematical Justification of a Thin Film Approximation 
for the Flow down an Inclined Plane 

Hiroki Ueno and Tatsuo Iguchi 
Abstract 

We consider a two-dimensional motion of a thin film flowing down an inclined plane under 
the influence of the gravity and the surface tension. In order to investigate the stability of 
such flow, we often use a thin film approximation, which is an approximation obtained by 
the perturbation expansion with respect to the aspect ratio of the film. The famous example 
of the approximate equations are the Burgers equation, Kuramoto-Sivashinsky equation, 
KdV-Burgers equation, KdV-Kuramoto-Sivashinsky equation, and so on. In this paper, we 
give a mathematically rigorous justification of a thin film approximation by establishing an 
error estimate between the solution of the Navier-Stokes equations and those of approximate 
equations. 


1 Introduction 


In this paper, we consider a two-dimensional motion of liquid film of a viscous and incompressible 
fluid flowing down an inclined plane under the influence of the gravity and the surface tension 
on the interface. The motion can be mathematically formulated as a free boundary problem for 
the incompressible Navier-Stokes equations. We assume that the domain 0(t) occupied by the 
liquid at time t>0, the liquid surface T{t), and the rigid plane S are of the forms 

( n{t) = {(x, y) G I 0 < y < /lo + r]{x, t)}, 

{ r(t) = {(x,y) G I y = ho+ y(x,t)}, 

[ S = {(x,y) G M^l y = 0}, 

where ho is the mean thickness of the liquid film and r]{x, t) is the amplitude of the liquid surface. 
Here we choose a coordinate system (x, y) so that x axis is pointed to the streamwise direction 
and y axis is normal to the plane. We consider fluctuations of the Nusselt flat film solution, 
which is the stationary laminar flow given by 

(1.1) yi = 0, Ml = (pysina/2//)(2/ioy - y^), mi = 0, pi = po - pg cos a{y - ho) , 


where y is a constant density of the liquid, g is the acceleration of the gravity, a is the angle of 
inclination, p is the shear viscosity coefficient, and po is an atmospheric pressure. Throughout 
this paper, we assume that the flow is Zo-periodic in the streamwise direction x. Rescaling the 
independent and dependent variables by using ho, lo, the typical amplitude of the liquid surface 
oQ) Uo = yy/i-o sina/2y, and Po = pghosma, the equations are written in the non-dimensional 
form 
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u = 0 on S, t > 0. 


Here, 5, e, R, and W are non-dimensional parameters defined by 
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where a is the surface tension coefficient. Note that 6 is the aspect ratio of the film, e repre¬ 
sents the magnitude of nonlinearity, R is the Reynolds number, and W is the Weber number. 
Moreover, we used notations u = {u,6v)^,u = {u,0)^,u = 2y — y‘^,Vs = {6dx,dy)'^, = 

Vs ■ Vs.Dsf = ^{Vsif'^) + {Vsif'^))'^}, and n = {-s6r]x,l)'^. In this scaling, the liquid 
domain and the liquid surface T^(t) are of the forms 


f ^eit) = {(a;, y)eR'^\0<y<l + ey{x, t)}, 
I ^eit) = {{x,y) G I y = l+er]{x,t)}. 


Concerning a mathematical analysis of the problem in the case of 5 = e = 1, Teramoto [14j 
showed that the initial value problem to the Navier-Stokes equations (ll.2h ~ (ll.4p has a unique 
solution globally in time under the assumptions that the Reynolds number and the initial data 
are sufficiently small. Nishida, Teramoto, and Win [10] showed the exponential stability of 
the Nusselt flat film solution under the assumptions that the angle of inclination is sufficiently 
small and x G T in addition to the assumptions in |14] . Furthermore, Decker |15] studied the 
asymptotic behavior for t —)• oo of the solution in the case of x G M and showed that the 
perturbations of the Nusselt flat film solution decay like the self-similar solution of the Burgers 
equation under the assumptions that the initial data are sufficiently small and R < Rc. Here, 
Rc = is the critical Reynolds number given by Benjamin [2]. On the other hand, Ueno, 

Shiraishi, and Iguchi [16] derived a uniform estimate for the solution of (jl.2|] - (jl.4|] with respect 
to 6 when the Reynolds number, the angle of inclination, and the initial data are sufficiently 
small. 

Benney [3| derived the following single nonlinear evolution equation 
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(1-5) r/t + 2(1 erjfr/x - tt(Rc - R)Shxx + CiS^rfxxx 
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C 2 e 6 (r]r]xx + vl) + ‘^-^Si^hxxxx = 0 ( 5 ^ + + £< 5 ^) 
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by using the method of perturbation expansion of the solution {u,v,p) with respect to <5 under 
the thin film regime (i <C 1. Here, Ci = C'i(R, a) and C 2 = 0*2 (R, a) are constants independent 
of 6, e, and W. Explicit forms of Ci and C 2 will be given in Section 3. Many approximate 
equations are obtained from (|1.5h by assuming that parameters e, W, and R have appropriate 
orders in 6. In the following, we assume £ = 5 and R < Rc and set 
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I. Burgers equation 

Assuming Wi <W <5“^W2 in (11.51) . we have 
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II. Burgers equation with a fourth order dissipation term 

Assuming W = 6 ~‘^W 2 in (|1.5I) . we have 

8 2 W 2 2 

Vt A rj^(Rc “t Sfjxxxx - 0(5 ). 
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Plugging (II.6p in the above equation and passing to the limit e = 5 —0, we obtain 

8 2 Wo 

(1.8) Cr + 4CCx - ^(Rc - R)Cxx + - Cxxxx = 0. 

15 3sma 

III. Burgers equation with dispersion and nonlinear terms 

Assuming Wi < W < W 2 in (11.51) . we have 


g 

rjt + 2r]x + ^eririx - —(Rc - R)5??xa; + CiS‘^r]xxx + C 2 e 6 {r]r]xx + vl) + ‘^'s'^v'^Vx = 0(5^). 

15 

Plugging (|1.6p in the above equation and neglecting the terms of 0(5^), we obtain 

(1.9) Cr + ‘^CCx ~ ~(Rc ~ R)Ca:a; + ^{Cl^xxx + C 2 {CCxx + Cx) 4" ^C^Ca:} = 0. 

IV. Burgers equation with fourth order dissipation, dispersion, and nonlinear terms 

Assuming W = 5~^W2 in (II.5|) . we have 

g 

r]t + 2r]x + 4.er]r]x - 7r(Rc - R)5r?xx 

15 

+ Ci6'^t]xxx + C2e6{7]rjxx + rjl) + 2£^rfr]x + ^-^^S'^rjxxxx = 0(5^). 

0 sin a 

Plugging (11.61) in the above equation and neglecting the terms of 0(5^), we obtain 

(1.10) Cr + 4CCx - 4 (Rc - R)Cxx 

15 
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We remark that (|1.9I) and (11.101) are higher order approximate equations to Burgers equation 
(II.7p . In this paper, we assume R <C Rc in order to use a uniform estimate in 5 for the solution 
of the Navier-Stokes equations, which is a severe restriction. Uniform estimates in 5 for the 
solution play a most important role in the justification for these approximation. Here if we 
could assume R > Rc, then (II.8p would be the Kuramoto-Sivashinsky equation (see |9], [12], 
and m)- If we could assume Rc — R = 5R > 0, then we would obtain the 5-independent 
KdV-Burgers equation 


(1.11) Cr + 4CCa; .. _ Cxx + CiCxxx — 0 

15 

by plugging ([OD in m and passing to the limit e = 5^ —>• 0 under the assumption Wi < 
W < W 2 . Moreover if we could assume Rc — R = —5R < 0, we would obtain the 5-independent 
KdV-Kuramoto-Sivashinsky equation 
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by plugging (|1.6I) in ()1.5I) and passing to the limit e = —)• 0 under the assumption W = 5 ^W 2 . 
More details or a list of useful references about the thin him approximation can be found in 

[g El El El [IB]. 

In this paper, we will give a mathematically rigorous justihcation of these thin him approx¬ 
imations by establishing an error estimate between the solution of the Navier-Stokes equations 
(ll.2p - (ll.4p and those of the approximate equations (I1.7l) - (ll.l0p . We note that we cannot just yet 
justify the Kuramoto-Sivashinsky equation, the 5-independent KdV-Burgers equation (jl.lip . 
and the KdV-Kuramoto-Sivashinsky equation (|1.12l) because without the assumption R <C Rc 
we have not yet obtain a uniform estimate in 6 for the solution. We also remark that Bresch 
and Noble [5] justihed the shallow water model by proving that remainder terms converges to 0 
as 5 ^ 0 (see also HD- 

The plan of this paper is as follows. In Section 2, we give our main theorem after we 
transform the problem in a time dependent domain to a problem in a time independent domain. 
In Section 3, we derive approximate solutions by using Benney’s method. In Section 4, we recall 
the energy estimate for the solution of the Navier-Stokes equations obtained in [16]. Finally, we 
give an error estimates in Section 5. 

Notation. We put 12 = T x (0,1) and F = T x {y = 1}, where T is the flat torus T = M/Z. For 
a Banach space X, we denote by || • ||x the norm in X. For 1 < p < oo, we put \\u\\lp = 
ll^ll = \u\lp = l)||/^p(']p), and |u|o = \u\i 2 . We denote by (•,-)n and (-j-jr the inner 

products of L^(n) and T^(F), respectively. For s > 0, we denote by and H^{T) the 

Sobolev spaces of order s on 12 and F, respectively. The norms of these spaces are denoted by 
II • ||s and I • Is. For a function u = u{x,y) on 12, a Fourier multiplier P{Dx) {D^ = —idx) is 
defined by {P{Dx)u){x,y) = T’(n)un(y)e^™*, where Un{y) = Jq u(x,y)e~^™^ dx is the 

Fourier coefficient in x. We put d~^f{x,y) = — f(x,z)dz and I2|/ = {(55x)*c^/ li + j = k}. 
f ^ g means that there exists a non-essential positive constant C such that f < Cg holds. 


2 Main results 


We rewrite the system (jl.2p -(jl.4p according to mm- Transforming the problem in the moving 
domain 12(2) to a problem in the fixed domain 12 by using an appropriate diffeomorphism, and 
introducing new unknown function [u',v\p') to keep the solenoidal condition, we obtain 
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where we dropped the prime sign in the notation and fi, f 2 , hi, and /12 are collections of nonlinear 
terms. See m for more details on the explicit form of these nonlinear terms. In the following, 
we will consider the initial value problem to (|2.1l) - (|2.3p under the initial conditions 

(2.4) n\t=o = m on r, (n,u)'^|t=o = (uo,uo)'^ in 0. 

Here, we assume r/o(x)dx = 0 and denote hi determined from initial data by h^^K 

We impose the following assumption on the non-dimensional parameters and initial data. 

Assumption 2.1. Let Rq, Ri, ceo) Wi, cq, and M be positive constants and m > 2 he an integer. 

(1) Conditions for parameters 

Parameters R, a,W, 5, and e satisfy 

Ri < R < Ro, 0 < a < oo, Wi < W, 0 < e = <5 < 1. 


(2) Smallness of initial data 

Initial data (r/o,uo,vo) and parameters W and 6 satisfy 

1(1 -h 6\D^\)‘^r]o\2 + 11(1 + |L'a;|)^('«o,'5uo)'^|| + 11(1 + \D^\)'^Dsiuo,dvo)'^\\ 

+ 11(1 + \D^\fDl{uo,6vof\\ + ,52w|( 1 + <5|Zl,|)r?oxl3 + VW||(1 + \D,\f6vo.y\\ < cq. 

(3) Regularity of initial data 
Initial data (r]o,uo,vo) satisfies 

11(1 + |.C)a;|)”*’''^(uo, -|- \rjo\m+4: < M. 

(4) Compatibility conditions 

Initial data {r]o,uo,vo) and parameters 6 and s satisfy 

{ uox + uoy = 0 in H, 

uoy + - 2(1 er]o)‘^r]o = ^ on T, 

rto = r’o = 0 on S. 

Remark 2.1. Under the assumption that there exist small positive constants Ro,ao, and cq 
such that Assumption 12 .1 1 is fulfilled, Ueno, Shiraishi, and Iguchi |16] proved the global in time 
uniform estimate with respect to 5 for the solution of the Navier-Stokes equations (|2.ip - (|2.4ll . 
See also Proposition 14.21 in this paper. 

For later use, we define the norm of a difference between the solution {rf^ ,u^,v^ ,p^) of the 
Navier-Stokes equations (I2.1I) - (I2.4I) and the solution C, of the approximate equations as 

(2.5) V{f,C,u,v,p) :=|7/^(t) - C(- - 2t,et)\l + ||(1 + mr{u^ - u)(t)f 

+ 11(1 + - t^)(t)f + 11(1 + \Dfir-\p^-p){t)\\\ 

where {u,v,p) is an approximate solution constructed from Q. Let ^ be the 

solution of (I1.7l) - (jl.l0p under the initial condition C|r=o = Vo, respectively. 

Now we are ready to state our main results in this paper. Note that the definitions of 
... appeared in the following statement will be given in Section 5. 
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Theorem 2.2. There exist small positive constants Rq and oq such that the following statement 
holds: Let m be an integer satisfying m > 2, 0 < Ri < Rq, 0 < Wi < W 2 , and 0 < a < uq. 
There exists small positive constant cq such that if the initial data (%, uq)'?^ o) the parameters 
6, e, R, and W satisfy Assumvtion \2. 11 then we have the following estimates. 

I. Burgers equation 

If the parameters 6 and W and the initial data tjq and uq satisfy 

(2.6) Wi<W<(5-1W2, \vo\m+7 + S-^\\{l + \D\^r+^uoyy\\<M CX:, 

then the following error estimate holds. 

(2.7) V{t-, v^,p^) < 

II. Burgers equation with a fourth order dissipation term 

If the parameters 6 and W and the initial data tjq and uq satisfy 

(2.8) W = r2W2, \m\m +12 + 6-^\\{l + \D\^r+\oyy\\<M <:x>, 

then the following error estimate holds. 

(2.9) V{t-,C^^,u^^,v^^,p^^) < C5^e-^K 

III. Burgers equation with dispersion and nonlinear terms 

If the parameters 6 and W and the initial data tjq and uq satisfy 

(2.10) Wi<W<W 2 , \r]o\m+i3 + 6-mi + \D\^r+\uoyy-ull^\t=o)\\<M < 00 , 
then the following error estimate holds. 

( 2 . 11 ) 

IV. Burgers equation with a fourth order dissipation, dispersion, and nonlinear terms 
If the parameters 6 and W and the initial data tjq and uq satisfy 

(2.12) W = r^W 2 , \vo\m+n + 6-^(1+ \D\,r+\uoyy-u^yl\t=o)\\<M < 00 , 

then the following error estimate holds. 

(2.13) V{t-< C5^e-^K 

Here, positive constants C and c depend on Ri,Wi,W 2 ,a, and M but are independent of S, s, 
R, and W. 

Remark 2.2. The assumptions for UQyy in (|2.6p and (|2.8I) represent the restriction on the initial 
profile of the velocity. Moreover, the assumptions for uoyy in (|2.10l) and (|2.12l) mean that the 
initial profile of the velocity have to be equal to that of the approximate solution up to 0(6“^). 

Remark 2.3. We see formally that the order of error terms in (11.71) is of 0{6), which implies 
that the error estimates (EH) and ()2.9p are natural. In a similar way, we see that the error 
estimates (|2.1ip and (|2.13p are natural. 

Remark 2.4. By introducing the slow time scale r = et, the norm decays exponentially and 
uniformly in r. 
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3 Approximate solutions 


In this section, following Benney’s perturbation method |3] we will give approximate equations by 
constructing approximate solutions. Hereafter, we assume e = 5. By straightforward calculation, 
we can rewrite (ISID-dSS]) as follows. 


(3.1) 

, 2 1 2 

6{ut + UUx + UyV) + —6px - —{5‘^Uxx + Uyy) = S — rjUyy + 5^/^ + 

xl xl xl 

6^(vt + UVx) + ^Py - ^S{6‘^VxX + Vyy) = ^‘^TlPy + 

^ '^X Vy - 0 


in Q, t > 0, 

in Q, t > 0, 
in H, t > 0, 


(3.2) 


6‘^Vx + Uy 

P- 5Vy - 


— 2(1 + 5?7)^r/ = 6^hi 
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- V H—^- "Hxx = 
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+ < 534 ^) 


on r, t > 0, 
on r, t > 0, 


(3.3) 


u = V = 0 on S, t > 0, 


(3.4) 


Pt + Px-v = (5^/i3 on r, t > 0 , 


where 


(3.5) 


= 1 {^P^Uyy - 2r]Px + 2yr]x;Py) + Ptu + yptUy 

+ y^PxU + 2y{y - l)pux -y‘^{y- 2)r]xUy - uux - vUy + 2{2y - l)r/u, 

/f ^ = 1 ( - 2p^py + 2yxUy + 2puxy) , 

( 2 ) 

h\- 2r}rjx; + PxU + pUx- 


We proceed to derive the approximate equations following Benney [3]. Let p = p{x,t) be 
a given function. For any 6 G (0,1], let {u,v,p) be the solution of (I3.1I1 - (|3.3I) and we expand 
(m, v,p) as 


(3.6) 


u = uq + 6ui + S‘^U 2 H-, 

< V = Vo + 6vi + 6‘^V2 H-, 

= Po + Spi + S'^Po H- 


and substitute these into (ISID-O, we obtain a sequence of perturbation equations for each 
order of 6. Here, uq and vq are different from initial data defined in (12.4p and hereafter we use 
this notation whenever it does not lead to confusion. By assuming W = 0(1), the 0(1), 0{5), 
and 0(5^) problems are as follows. 


^Oyy — 0, 

< UOy = 2p, 


POy — 0, UQx 4“ VQy — 0 
1 

Po = - - P 

tan a 


uo = vo = 0 


(3.7) 
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in Q, 
on r, 
on S, 









(3.8) 


(3.9) 


Uiyy = R(uot + (2y - y^)uQ^ + 2(1 - y)vQ) + 2pox + 27/7/0^ ™ 0, 

2piy = VQyy + 2r]P0y, Ulx + Vly = 0 111 0, 

Uiy = Aif, Pi = -uqx on r, 

tti = vi = 0 on S, 


1 ( 

U2yy = R(wn + (2y - y^)uix + 2(1 - y)vi) + 2pi:c + 2r]Uiyy - uq^x - R/i ’{r],UQ,vo,po) 

2.p2y = Vlyy + 2ripiy - R(7;ot + (2y - y‘^)vox) + R/Piy, uo,vo,po), U2x + f2j/ = 0 
U2y = -Vox + 2r]^, P2 = -Uix + {y, Uq) - 


W 


sma 


~^xx 


U2 = V2 = 0 

Solving the above boundary value problem for the ordinary differential equations, we have 


in n, 
in n, 

on r, 

on S. 


(3.10) 


no = 2yy, 

Vo = -y‘^yx, 

1 

po = I — y, 

tan a 


(3.11) 


'ui = {^y^ - y) Ryt + {(y^ - 2?/)^ + (i/ - fy) R} yx + 

= {-T2y^ + hy‘^)^p^t + {{-h^ + y^)tic^ + {-my^ + h‘^)^}p=ox-^y‘^yyx, 


^Pi = -{l + y)yx, 

(3.12) 

U 2 = {^y^ - + ^y) R^ytt 

+{(^2/^ - + ly)i^ + - ly^ + 


+{ {-ly^-y^ + 5y) + {-±y^ + ±y^-ly^ + ly) 


R 

tan a 


+ (-5bo2/® + 315^^ - ^y^ + is?/) 

+2yy^ + R(|y3 _ 4:y)yyt + {R(/ - Ay) + {3y^ - Qy)t^}yyx, 

^2 = {-my^ + iiy^ - ^y^) 

+ { i-my^ + T2y^ - W) + i^y^ - ai?/^ + my^ + - my^) 

+{(1/ + - hi + (io?/^ - - bl 1^ 

+ (sio?/^ “ TsW?/* + ^y^ ~ T^yl ^‘^}yxxx 

-^yWyx + R( - + 2y‘^){yxyt + yytx) + {R( - b^ + 2y^) + {-y^ + ^yltb^}(yx + yyxx), 

P2 = ih + s) ^y^t + + bb^ + y + h)tbA + bwy^ + b^ + b + w) b y^^ 

+ {R(4y - 4) - 5y + 3}yyx. 

Using the above expressions, we put 

T T T / \ T T T / \ 

= Po, 

= P1, 

= P2- 

In view of the perturbation expansion (13.61) substituting v = + 6v{^^ + into (13.4p . we 

obtain the approximate equation 




'uo^{y;y) 

= Uo, 

Vo^ly;y) 

= Vo, 

Poly,y) ■ 

(3.13) 

< 

ui^ly;y) 

= Ul, 

vi^ly;y) 

= Vl, 

Pily,y) : 



,«P(2/;h) 

= U2, 

V2ly;y) 

= V2, 

P2ly,y) : 


yt + 2yx + Aeyyx “ Y^(^c - R)Syxx + Ci5 yxxx + C2e5{yyxx + yt) + 2e y yx = 0{5 ), 






















where Ci = 2 + IR^ _ 40 R 


and 6*2 = — 


63 63 tan a ^ 5 tana* 

Thus far we have assumed W = 0(1). Taking into account that W is contained only in the 
second equation in (j3.2p and modifying 0{6) problem under the assumption W < 0(<5“^), we 
see that {uq,Vq,Pq) and which are defined by 


(3.14) 


f «o(y; V) ■= uq, vl{y; tj) := vq, pl{y, p) ;= po, 

\u{{y,p) := ui, v{{y;p):=vi, p{iy,v) ■= Pi --^Vxx, 

are the solutions of the problems. Putting v = Vq + dv( and substituting this into (13.41) . we 
obtain the approximate equation 

g 

pt + + Aeppx - —(Rc - R)Spxx = 0{6'^). 

15 

Similarly, modifying 0(1) and 0(5) problems under the assumption W = 0(5“^) and putting 

13.151 ho-=^ 

■■= ui - ^{y^ - 2y)p, 
we obtain the approximate equation 


:= Vo, Po := Po 

5'^w I 


'-^Pxx, 
siriQ 


uf := vi + - y‘^)pxxxx, Pi ■= Pi, 


8 2 W 2 

Pt T ‘^Px T ^^PPx Ty (Rc R)5r/a;a: T ^Pxxxx 

15 3sma 


= 0(52). 


Moreover, putting 
(3.16) 


..iv 

Uq 

ui^ 


= Uo, 

= Ul, 

= U2- 


IV 

Vo := Vo, 
IV 

Vi := vi, 

m 


IV 

Po ■=Po, 

IV 

Pi ■= Pi - 


<5W 


Pxx, 


(7/2 - 2y)px 


JV _ 


:= V2 + 


5VI /1„.3 


{\y^ - y^)pxxxx, pi^ ■= P2 + 


w 


Vxx 


and V = Vq -\-6vI +5‘^V2 and substituting this into (j3.4l) . we obtain the approximate equation 


Pt + 2px + Aeppx - — (Rc - R)Spx 
15 


2 W 2 


T Oi5 Pxxx T C2£5{tjPxx T Px} T 2£ p Px T „ . 5 Pxxxx — 0(5 ) 


3 sin a 


under the assumption W = 0(5 


-l^ 


4 Energy estimate 


In this section, we will derive energy estimates. Let p = p{x, t) be a fixed function and (u, v,p) = 
{u{y, p),v{y; p),p{y; p)) be an approximate solution constructed from p satisfying Ux+Vy = 0 and 
u\y=o = v\y=o = 0) which will be defined in the next section. Using the approximate solution, 
we define ilii,'il) 2 , 4 >i: 4 ’ 2 , 4>3 by the following equality. 

ipiiy^p) := + vux +Uyv) + ‘^^Px - ^{S'^Uxx + Uyy) - 5f\^\p,u,v,p)^, 

'^p2{y,p) ■= ^^5‘^{vt + UVx) + “^Py - ^S{5'^Vxx + Vyy) - 6f^^\p,U,p)'^, 

+ '^y~ 2(1 + ^'nf'n]\y=i-, 


(4.1) 


53 

1 




Mp) ■= pivt + Px-v- 5‘^ho{p)}\y=i, 


—p + A^Pxx - I 

n a sm a } 


y=l 
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where 


( 4 . 2 ) /f ^ = -^VUyy + - 5 /f ^■ 

Then the approximate solution satisfies the following equations. 

(4.3) 

^ 2 1 

6{ut+ uux+ Uyv) +—6px -—{6‘^Uxx+ Uyy) = Sf[^\r],u,v,p) + S^-ipiiy;!]) in Q, t > 0, 

2 1 

6‘^{vt+ uvx) +—Py -—6{6‘^Vxx +Vyy) = 6f^^\r],u,p) + 6^'4!2{y;p) in Q,t>0, 

rl sx 

^ Ux + Vy = 0 in n, t > 0, 


(4.4) 


6 ‘^Vx + Uy — 2(1 + 6 p)‘^p = 6 ^(l)i{rj) on T, t > 0, 

' p - 6 vy -+ ^r^rjxx = S‘^hf \r],u) + 6 ^(j) 2 {p) on T, t > 0, 
tan a sin a 

. Vt + Vx-v = 6 ‘^h 3 {p) + 6 ^(l) 3 {p) on T, t > 0, 


(4.5) 


u = V = 0 on S, t > 0. 


Note that in the next section, we will give explicit forms of ipi, ip 2 , 02; and 03. Let 
{p^,u^,v^,p^) be the solution of (I3.1I) - (I3.3I) and we set 

H := p^ — p, U := — u, V := — v, P := p^ — p. 


Taking the difference between ()3.ip ~ ()3.4p and (|4.311 ^ (14.5p . we have 


(4.6) 


(4.7) 


6{Ut + uUx + UyV) + '^^P'x ~ '^{^‘^Uxx + Uyy) 

= Fi + 6^f[^\p^,u^,v^,p^) — 6^'il>i{y;p) in fl, t > 0, 

' 6 HVt + uVx) + lPy-^ 5 { 6 ^Vxx + Vyy) 
rl sx 


= F 2 + 5^f2^\p^,u^,v^,p^) - 5^'ip2{y]p) in fl, t > 0 , 

Ux + Vy = 0 in n, t > 0, 

6‘^Vx + Uy — (2 + b{p^, p))H = 6^hi{p^,u^,v^) — S^cj)i{p) on T, t > 0, 

' P - 5Vy - ^ — H + \ ^ Hxx = G 2 + 5^h'"2\p^,u^,v^) - 6^(j)2ip) on T, t > 0, 

tan a sm a 

Ht + Hx-V = G^-5^3{r]) on T,t> 0, 


(4.8) 


U = V = 0 on S, t > 0, 


where 


(4.9) 


'Fi = 6{f[^\p\u^,v\p^) - f[^\p,u,v,p)), F 2 = 6{f^^\p^,u\p^) - f^^\p,u,p)), 
< b = 25[6{p^)‘^ + (2 + 6 p)p^ + Sp"^ + 2 p ), 

G2 = 6‘^{h^2\v^,u\v^) - h‘'2\r],u,v)), G3 = 5 '^{h 3 {p^) - h 3 {p)). 
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For convenience, we set 


U:={U,6Vf, F:={F,,F2f, /:=(/P,/f)T, r/; := (V^i, 

We proceed to derive an energy estimate to (j4.6j) - (j4.8[) following [16]. In view of the energies 
obtained in [T6| (see (3.6)~(3.8) and (3.24) in [l6]), we put 


MH,U) 




tan a 


sin a 


^x\o 




Sin a 


^xx\{) 


[ R \ tan a sin a J ] 

+ fhis^utf + 4( —< 5 = 111 , i;? + '= 


R \ tan a ” sin a 

^o{H, [/, P) := 5\\U^f + 5\\d-^P^f + 6\H^\l + + 5^W^\H, 

+ 6\\VsU,f + 6^\\VsU^A^ + 6\\VsUtf 


2 

XXX 10 


Here, / 3 i,/? 2 , and /Js are appropriate positive constants (see (3.28) in [E]). Integrating by parts 
and using the third equation in (14.711 and Poincare’s inequality, we see that for any e > 0 there 
exists a positive constant such that 

5^\{{F + 6^f - U,,h\ < + CMF.f + <5®ll/xf + 

\{H, {bH),)r\ < eS\H,\l + C.b-^bHUl 

(677)x)r| < eS^W\H^^\l + C,6W\{bHUl 

5^W\{H,„ Gs - (^303)r| < e5^W\H^^\l + C.i^WdGslg + (^®|<^3lo), 

(^®W|(i7,,,,,d3<^3xx)r| < e<5^w2|i7,,,|2 + C,(5i3|^3xxx|o, 

S^W\{H^^t,G3t - 6^ht)r\ < + 

+ Q (1 + W2)(i3(|G3dg + 6^\<P3t\l) + SH\G3a^x\l + S^\^3xx\l). 

Here, we used the inequality |R(-,l)|o = |R(-,1) — R(-,0)|o < WVyW = \\Ux\\ thanks to the third 
equation in (|4.6I) and the second equation in (|4.8|) . In the following, we use frequently this type 
of inequality without any comment. Taking into account the above inequality and (3.27) in jl 6 ] . 
we need to estimate the following quantities. 

(4.10) PPoHZ,) := (5W + 5-^)\ibH)x\l + S^\{bH)xx\l + 

+ (^ ^|G2|o + <5|G2x|o + b‘^\\Px\^G2x\o + <^1(^21) I 
+ 5W|G3|g + 5"|G3x|g + S^lGsxxll + 5"w2|G3i|g + 5^W\{Hxxxx, G^xxhl 
+ 6-^\\Ff + 6\\Fxf + 6\{Ft,Uth\, 

(4.11) pZo^{Z 2) :=S^\hi\l + 6'^\hix\l + S^\\Dx\'^hix\l + S\hu,Ut)^ 

+ S^ff + 6^\\fxf + S^\{ft,Uth\, 

(4.12) =ylQ^(Z3) := 5 ®|(^i|q + 5'^|(^12;|q + 5®||Zlx|2(()ix|o + <^'^|'?^iilo + <5^|02|o + '^^l'/’ 2 x|o 

+ 11 .Da; I 2 (j() 23, I g + (5'^|(()2 i|o + 5'^W|i;/)3|g + 5^|((>3x|o + b^^\4^3xx\o 

+ ^^^\4>3xxx\o + ^^W^|(/>3do P + '^'^IIV’xll^ + 
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where 


Z, = iH,U,bH,G 2 ,G 3 ,F), Z 2 = {U,h,h^i\h3,f), Z 3 = <^ 2 ,03, V’)- 

For an integer m > 2, we set 

m m 

(4.13) >^m{H,U) ^m{H,U,P) ■.= Y,-^o{d^H,d^U,d^,P), 

k =0 k =0 

m 

(4.14) ^J,{H,U,P-,v) := + 

k =0 

m 

(4.15) pC{U) :=Y.^A\d^^Z 2 ), 

k =0 

m 

(4.16) ■.= Y,{^o\dlZ^) + \{dlH,5^dlcl2:,)r\]. 

k =0 

Here, the terms Yl'k=o ^xG' 3 )r| and Yllk=Q K^x^,'^^^x'^ 3 )r| come from (3.30) in [T6j. Ap¬ 

plying to (|4.6p - (|4.8p . using [16, Proposition 3.2], and adding the resulting inequalities for 
0 < k < m, we obtain the following lemma. 

Lemma 4.1. There exist small positive constants Rq and ao such that i/0 < Ri < R < Rq, 
Wi < W, and 0 < a < oq, then the solution (H, U, V, P) of (|4.6l) " (|4.8p satisfies 

(4.17) + + 

where the constant G is independent of 5, R, and W. 

For later use, we modify the energy and the dissipation functions (om and as 

(4.18) 4n(77, U) := ^m{H, U) + ||(1 + \Dfiruf + ||(1 + \D^\rUy\\\ 

(4.19) ^m{H, U, P) := U, P) + <5|(1 + + i 6 ‘^WfS^\\Dfi" 2 H\l, 

+ <5-^(1 + |71.|)”*(1 + <5|Z1^|)(V5P, Uyy)f + ^11(1 + . 

We also introduce another energy function by 

(4.20) 

U) :=!(! + + 52||(i + |I1,|)-R||2 + 52||(i + \D^\)--U^f 

+ 11(1 + \D^\rDluf + (<52w)2|(1 + ,5] 11,1)77,1^+1 + VW||(1 + \D^\r5V^y\\\ 

which does not include any time derivatives. Setting Em = ^m{v^,u^) and Fm = ^mih^,P^) 
and using [16, Theorem 2.2 and Proposition 6.1], the following uniform estimate holds. 

Proposition 4.2. There exist small positive constants Rq and oq such that the following state¬ 
ment holds: Let m he an integer satisfying m > 2, 0 < Ri < Rq, 0 < Wi < W 2 , and 0 < a < oq- 
There exists small positive constant cq such that if the initial data {r]o,uo,vo) and the parame¬ 
ters 6 , e, R, and W satisfy Assvmption \2. 1\ and W < 6 ~^'W 2 , then the solution {r]^,u^,v^,p^) 
of (j2.ip - (|2.4p satisfies 

roo 

E 2 {t)<co, supEm+i{t)+ Fm+litW < G, Em+i{t) < 
t>o Jo 

Here, positive constants G and c depend on Ri,Wi,W 2 ,a, and M but are independent of S, e, 
R, and W. 
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Moreover, we easily obtain the following lemma. 


Lemma 4.3. Let a > 0, 0 < Ri < R < Re- There exists small positive constant ci such that if 
s > 2 and \ri 0\2 < ci, then the problems ()1.7ll -- (ll.l0p under the initial condition C|t=o = Vo have 
unique solutions , and , respectively, which satisfy 

poo 

sup|C^(r)|2+ / |Ci(T)|2dr < Clryol^, |C^(t)| 2 < C'|ryo|^e"^'^‘, 

T>0 Jo 

poo 

sup |C^^(r)|f + / {\Ci^{T)\l + lCi^(T)|2)dr < C\r]o\l, \C^\t)\1 < C'|r/o|^e“'^'^*, 

T>0 Jo 

poo 

sup|C^^^(r)l2 + / |C"^(r)|fdT < Clr/ol^, < C|r/o|«e■'=^^ 

T>0 Jo 

poo 

sup|C^^(t)| 2 + / (|Ci^(r)|2+ 5|Cir(^)ls)dr < C|r?olL |C-^^(t)| 2 < C|r/o|^e"'=^*. 

T>0 JO 

Here, Rc = is the critical Reynolds number and positive constants C and c are independent 

of 6 and R. 


5 Error estimate 


We will show ()2.1ip under Assumption [24] and p2.10l) . We can show the other claims in Theorem 
12.21 in the same way as the proof of (I2.11h and we will comment about the difference at the end 
of this section. Let be the solution of (II.9p under the initial condition C^^^\r=o = Vo 
we put ri^^\x,t) := — 2t,et) and 




( III f 
u (x, 

y,t) 

(5.1) 

< 

v^^^{x, 

y,t) 



y^^{x, 

y,t) 

where 

yin 

“0 

III in 

Oq ,Po ,■ ■ 

(5.2) 


vi^^ 



(x, y, t) := r]^^^{x, t)) + 6u{^^{y, t]^^^{x, t)) + 5‘^u{^^{y, r]^^\x, t)), 

v^^^fx, y, t) := ul^^{y; rj^^^x, t)) + 6v{^^{y; rj^^^x, f)) + d^v^^^{y; rj^^^x, f)), 
jn„in+\\ I a^iii(... ^111 1 x 2 ^iii^111 




+ —(Rc-R)<5r?3,3, -Ci6p. 


2^111 

XXX 




Using the approximate solutions p5.ip . we define 'ifi, 1 ^ 2 , 4‘i, 4>2, and cf^ by dn]). By using the 
equality (15.2p to eliminate the t derivatives of v^^\ can rewrite these terms as follows. 


(5.3) 


''ifi{y,p^^0 = Ciiy)dW^^ + C2{y)ddfr]^^^ H- ^ C7{y)6^d^r]^^^ + , 

V’2(y; V^^^) = C?,{y)dW^^ + C9{y)6d^p^^^ H-h Ci5{y)6^d^^r]^^^ + 

< + • • • + C2id^dy^^ + Ni^^, 

= C22dW^^ + C23SdWii + • • • + C2ed^dW^^ + 

= C27dW^^ + C235dW^^ + • • • + Csod^dW^' + 


where Ci,... ,Ci 5 are polynomials in y, Cie,... ,^30 are constants, and ,..., are col¬ 
lections of the nonlinear terms of the form 


(5.4) 


1 
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Here we generally denote polynomials of f by the same symbol = ^{f) and $0 is such a 
function satisfying <l>o(0) = 0. We also use such a function $0 depending also on y E [0,1] and 
denote it by <ho(/;y)) that is, <ho(0;y) = 0. Let {rj^,u^,v^,p^) be the solution of (I‘2.1|) - (I2.3|) and 
we set := rj^ — := (u^ — , 6 {v^ — := and so on. 

We prepare several lemmas to proceed the error estimate. 

Lemma 5.1. Under the same assumption as Proposition \4-^ for any e > 0 there exists a 
positive constant such that we have 

(5.5) jV^{U^^\t) < e^mit) + C,5^Em{t)Fm+l{t), 
where is the collection of nonlinear terms defined by ()4.15l) . 

Proof. By the explicit form of f, hi, and (see (|3.5p and Section 2), we can obtain the 
desired estimate in the same but more easier way as proving [16, Lemmas 5.11 and 5.12]. □ 

Lemma 5.2. Under the same assumption as Proposition 14-31 for any e > 0 there exists a 
positive constant Ce such that we have 

where is the collection of nonlinear terms defined by ()4.16l) . 

Proof. By the well-known inequalities ||5^(/y)|| < ||/||i,°° ||5^5ll + lbl|L°° ll'^x/ll and \\d^^o{f-,y)\\ < 
C{\\f\\L--)\\d^f\\, dOIl-dSai) leadtoX]r=o^^(^x^3) < {^ + W^^\‘L+i 2 )^^\hx^\‘L+i 2 - Moreover, 
by Poincare’s inequality and (15.4p . we see that \{d^H, 6 ^d^ 4 > 3 )r\ < e 6 \d^Hx\l + < 

+ Cf:{l -|- W^^\‘in+i 2 )^^\'nx^^\‘in+i 2 - These together with Lemma [TS] imply the desired in¬ 
equality. □ 

Lemma 5.3. Under the same assumption as Proposition 14-31 for any e > 0 there exists a 
positive constant such that we have 

(5.6) 

<iC\E 2 {t) + e)#"^(t) + C,{Emit)^i'^{t) + 5^Emit)Fm+i{t) 

+ i't)\m+l2 + {Fm{t) + ^Wx\t)\^+I2)^m^{^)]-: 

where is the collection of nonlinear terms defined by ()4.14l) . 

Proof. In this proof, we omit the symbol III appeared in a superscript of solutions for sim¬ 
plicity. By (13.5p . (14.21) . and (|4.9I) . we see that F is consist of terms of the form 

' 6 Mv\Svi-,y)i^SUy,VsP) + 6 \v')HUyy,Py), 

6Mv\s4,u^-,y)i6V,dUx), 

< 6 <I>o{Sy^^,dvf,dv^-y)(U,Uy), 

6 ^o{p, u, Vsu, Vsuy, Vsp; y){dH^, 6 Ht, U, 6 V), 

+Py)H 

and that G 2 = S‘^{r]^{ 2 Hx + Ua;) + piU -h { 2 r]^ + Ux)H PuH^}, G 3 = S‘^{{r]^)‘^Hx ++ p)rjxH}, 
and bH = 2S[S{p^)‘^ -|- (2 -|- 6 p)p^ -|- -|- 2p)H. Note that using (15. 1|) and (15.2p , we can express 


14 









the approximate solutions u,'Vsu,Uyy, and Vsp in terms of r] and its x derivatives. In view of 
these, by putting 


$1 = ^(r]\ Spl, 6r]f, v), 

$4 ^ ^{y^5yx,---,S^^dl°y;y), 


'W := (SHx, 6Ht, 5‘^Hxx, S^Htx, 6V, dUx,SUt, 5VsUx, dVsUt, VsUy, VgUxy, 

< VsP, VsPx, SUx\r, SUtIr, S^Uxx\r, 6^/‘^\Dxf/^U\r), 

:= {H,6Hx,6Ht,S^Hxx,6‘^Ht^,S^Hxxx,U,VsU,6Ut,U\r), 

it suffices to estimate 

'h = 5\\d^,{^lW)r, 
i2 = d\\d^,{^lQ)r, 

h = S^\{d^{v^Utx),d^Vt)r\, 

^ h = 6^\id^,i‘^lVsUty),d^Uth\, 

' h = 6^\id^,i<^lVsPt)d^xUth\, 

h = S^^\id^Hxxxx,d^G3xx)r\ 

for 0 < k < m. 

By Proposition 14.21 and ||(m,u)||l°o ^ IK^yi'^j/)ll + \\{uxy,Vxy)\\ thanks to the boundary con¬ 
dition u\y=o = v\y=o = 0, we obtain 

(5.7) \\H\\U<E2, \\d^xHf + \\d^x^lyf<Em, 

(5.8) \\^l\\l^<F2, \\d^x^lf + \\d^^^lyf<Fm. 

In the same way as the proof of Lemma 15.21 we have 


(5.9) 6\\^ 


< 


sh 




md^xKi 


2 


0y\ 


< 




x\m+l2J 




4|2 

Olm— 



2 

m-|-12. 


On the other hand, it is easy to see that 


(5.10) IllPf + IllP.f < #2, Wd^xWf < #, 

( 5 . 11 ) \\Qf + \\Qxf<^2, \\d^xQf<<^m, 


where we used the trace theorem |/|q -|- ( 5 ||L) 2 ;| 2/|2 < ||J||2 _|_ (5^||/3,||2 -)_ ||/y|p to estimate the 
term 2 [/|q. In the following, we often use the inequality 

(5.12) ||8j(a/)|| < ||a||i»||8j/|| + (||8ja|| + ||8Xll)(ll/ll + IIMI). 


which have been shown in [16, Lemma 5.2]. 

As for Ii, by ()5.7I) . (|5.10p . and (15.121) . we have Ii < E 2 ^m + Em’^ 2 - As for I 2 , by (|5.8I) . 
(|5.1ip . and (I5.12h . we have I 2 < Fm^m- As for I3, by integration by parts, we have < 
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CeS^\r]^Utx\l^_i + 1 < C'e(^ 2 ^m + Em^ 2 ) + As for I 4 , by integration by parts 

fn 2 2 

in y, we have 

h < CJ^{\m^lVsUt)f + \m^lyVsUt)f) 

< -^ 4,1 + -^ 4,2 + -^ 4,3 + e^m.) 

where /4,i = a< 52 (|| 5 ^(ch 3 V 5 [/i)f + ||9^(cJ>3^V5[/i)||2),/4,2 = S^\{d^{<^lUtx), d^Ut)r\, and 14,3 = 
6 ‘^\{d^{^QUty), d^Ut)r\- The estimates for 14^1 and 14^2 are reduced to the estimates for Ii and 
Is, respectively. Thus, taking into account that we can eliminate the term ?7y|r in 14^3 by 
the first equation in (14.7p . this together with the estimates for I 2 , I 3 , 6 ^hi, and yields 
I 4 < e^rri + Cf:{E 2 ^m + .^m(^2 + + \'n\m+ 12 ^^\Vx \m+i 2 )}- As for /s, it suffices to show 

the case of A; > 1 because we can treat easily the case of /c = 0. Integrating by parts in x, 
and dSH]), we have h < e 6 ^\\d^Utxf+ CMd^~H^o^5PtW < + ^^^ 2 )• As 

for Ig, by (j5.7p . (j5.9|] . (j5.1ip . and (I5.12|] . we have 

h < 6{\\^l\\l^{\\d^,<^Y + )(IIQf + WQxf) 

+ )(||Qf + lig.f) + ll‘i>'llioc||<h4||i^||5,"gf} 

^ {Em + \'n\m+ 12 )^\Vx\m+ 12 ^rn- 


As for Ij, it suffices to show the case of k > 1 because we can treat easily the case of A; = 0. By 
the third equation in (14.71) . integration by parts, and the trace theorem, we have 


I 7 < CJ^\\Dx\-^dr\<^t>yt)\o + + CJ^\d^dMo 

+ e{5^\\Dx\"^d^,Vt\l + 6^\d^,Vt\l) 

< It,! + It ,2 + It ,3 + 


where/7,i = 17,2 = CJ^\d^{<f>f,Hxt+^^oG3t)\l and 17,3 = 

By the trace theorem, the second equation in (14.6p . and p5.9p . we have 


C,6^\5^dMl 


l7,i < 

< \<\l_iS^\\Utxxf + 5\^l\l^i5^\\d^xUtf + S^\\d^,Vtf) 

^ \v\m+12’^2 + S\7lx\m+12<^m- 


Recalling the explicit form of G 3 , we see that the estimate of It ^2 is reduced to Iq. Taking 
into account that we have already estimated 17^3 in the proof of Lemma 15.21 we obtain Ij < 
Ce{\'n\m+ 12’^2 + {Em + Hm+ 12 ¥\Vx\m+ 12 ^m + |?/lm+12<5^|m +12 } + • As for h, integration 

by parts, (IS21), cinci (j5.9p Iccid to 

5^W\{d’:,Hxxxx,d’^,G3xx)v\ < e{5^Wf5^\\Dx\'^H\l + G,6^\\Dx\^G3\l 

< + <^|lm+ 12)^2 + E2^m} ■ 


Therefore, by the boundedness of the terms Em and |??|m+i 2 which comes from Proposition 
M\ and Lemma 14.31 the proof is complete. □ 
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Lemma 5.4. Under the same assumption as Proposition \4-^ we have 


(5.13) 

+ 5^{Em+i{t) + |r/^^^(t)|™+i2), 

(5.14) 

^ ^(^) + {Fm{t) + 5\'nx\t)\l^+l2)^^\t) 


+ 6'^Em{t)Fm+l{t) + d^\'nx\t)\l^+l2, 

(5.15) 



Proof. In view of the discrepancy of non-homogeneous terms in the equations, modifying the 
proof of (6.2) in [16, Lemma 6.2], we obtain (15.1311 . Taking into account that we can eliminate 
Uyy in by using the first equation in (14.61) . modifying the proof of (6.3) in [16, Lemma 6.2], 
it is not difficult to check that (15.141) holds. Moreover, modifying the proof of (6.10) in [TB], we 
obtain (|5.15p . □ 


Lemma 5.5. 


Under the same assumption as Proposition ^.^ we have 



( 0 ) < 


Remark 5.1. This lemma together with ()5.15p yields 


(5.16) 


4l''(0) < 


Proof. By the second and third equations in the compatibility conditions, we see that 


(5.17) uo{x,y) =yuoy{x,l) - 


y /■! 


0 Jz 


UOyyix, w)dwdz 


{2yrio + 4y(5?7o + +dy{- 6vox + ^) - [ [ uoyy{x, w)dwdz 

Jo Jz 


It follows from (I2.10p and ||(1 + \Dx\)^^^Uyy^\t=o\\ < 6 (see the explicit form of that is, 
(I3.10p - ()3.13p and (|5.ip l that ||(1 + iDxD'^^moyyW < 6 . Thus, by (|5.17l) . the explicit form of 
()2.10l) . and the uniform estimate for 6 ‘^\hf^'^\m+i (see the proof of Lemma l5.II) . we obtain 
II(1 + \Dx\)'^^^U\t=o\\ < 5. Combining this and the first equation in the compatibility conditions 
leads to ||(1 + |T)j,|)”^17|j=o|| ^ Therefore, in view of the definition of (see (14.201) 1. using 
these and 77|j=o = 0, we obtain the desired estimate. □ 


Proof of (12.lip in Theorem \2.2[ By Proposition 021 Lemmas ITTl 15.11101 and (|5.13l) and (I5.14p 
in Lemma 15.41 if cq and e are sufficiently small, then we have 

(5.18) < Ci{ipi{t)S’l^^{t) + Em{t).^i^\t) + 5'^^p2{t)), 
where 

(5.19) (pi{t) = Fm{t) + <5|r?"^(t) 1^+12, Mt) = ^m(t)^m+l(t) + S\Vx^{t)\‘L+12- 

By considering the case of m = 2 in (|5.18l) and using Gronwall’s inequality and Proposition 14.21 
if Co is sufficiently small, then we have <^ 2 ^\t) + Jq ^ 2 ^^(s)ds < ^psit), where 

(5.20) (fsit) = S' 2 ^\ 0 )ex.p (^Ci J ipi{s)ds'^ + Ci J (5"^(^2('S) exp i: (^l(cr)dcr^ ds. 
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which leads to 


(5.21) f ^^^^{s)ds < ip3{t). 

Jo 

Note that by Proposition 02] and Lemma r4.8[ we have the exponential decay estimate for Em+i{t) 
and This together with (|5.18p . Gronwall’s inequality, and which 

comes from \H\o < \Hx\o and ||P|| < \\Vy\\ = ||?7a;|| (see (I4.13h and (I4.19p l yields 

< |(f^^^(0)exp (Ci ¥?i(s)ds^ + 

where 

(5.22) (p4{t) = Cl j (^2^^('5) + <^^Tm+i(s)) exp ^Ci / (/9i(cj)dcr^ ds. 

Combining the above inequality and p5.13p and (I5.15P in Lemma 15.41 we obtain 

(5.23) < C2{5^ + <^^(0) + 

Here, recalling the definition — 2t,et) and the assumption e = 5 and using 

Lemma 0l3l we have /o°° (^l^x^^(t)lsdt = ^ /o°° (^lCi^^('^)lsdr < \r]o\s. By this, the integrability of 
Fm+i which comes from Proposition 14.21 and p5.16p . we have ^psit) < 5^ (see p5.19p and p5.20p L 
This together with (I5.2ip leads to (p 4 :{t) < 5^ (see (|5.22p l. Combining this, p5.23p . and Lemma 
[531 we have 

(5.24) < Cs^^e—^ 

which implies (see (|2.5p and (14.181) 1. Here, we used ||H|| < 

||Py|| = ||Ha;||. Moreover, by taking into account the equality P(x, y, t) = P{x,l,t)—Jy Py{x, z,t)dz 
and using the second equation in (|4.6p . the second equation in (|4.7D . and the uniform esti¬ 
mate P5.24I) . we easily obtain ||(1 -|- \Dx\)'^{p^ — p^^^){t)\\‘^ < Note that in the case 

of 0((5“^) < W < 0{6~^) we can estimate the term which comes from the second 

equation in (|4.7p by <^^+ 1 - Therefore, the proof of p2.1ip in Theorem 12.21 is complete. □ 

We proceed to prove (12.71) . (12.91) . and (I2.13p . Let , and be the solution for (|1.7I) . 

pi.8h . and (ll.lOh . respectively under the initial condition C^|t=o = C'^^|t=o = C^^\t=o = Vo- We 
put rj^{x,t) := C^{x — 2t,et), r]^^{x,t) := C,^^{x — 2t,et), r]^^{x,t) := {x — 2t,et) and 

u^{x,y,t) := u(){y,r]\x,t)) + 6u{{y;r]^{x,t)), 

< v^{x,y,t) := u(^{y;p^{x,t)) + 6v{{y,y^{x,t)), 
y{x,y,t) := p(,{y;y^{x,t)) + 6p{{y;p^{x,t)), 

u^^{x, y, t) := u()^{y; y^\x, t)) + 6u{\y; y^\x, t)), 

< v^^{x,y,t) := {y;y^^{x,t)) + 5vi^{y,y^^{x,t)), 

y^{x,y,t) :=pl^{y;y^^{x,t)) + 6p{^{y;y^^{x,t)), 

'u^^{x,y,t) := n^^(y;7?^^(x,t)) -h5u{^(y;yf^(x,t)) -h f)), 

< v^^{x,y,t) := u()^{y,y^^{x,t)) + 5v{^{y;y^^{x,t)) + {y;y^^{x,t)), 

y^{x, y, t) := pi^{y, r?^^(x, t)) + 6p{^{y, r/^^(x, t)) + 6‘^pi^{y, r/^^(x, t)), 

where Uq,Vq,Pq, ... were defined by p3.14p ~ (l3.16p . In view of this, by applying the same argument 
as showing (12.111) . it is not difficult to check that (12.71) . p2.9p . and (I2.13P holds. Therefore, the 
proof of Theorem 12.21 is complete. □ 
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